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ABSTRACT 

We present multicolour images of the hosts of three z — 2 QSOs previously detected 
in i?-band by our group. The luminosities, colours and sizes of the hosts overlap with 
those of actively star-forming galaxies in the nearby Universe. Surface brightness radial 
profiles over the outer resolved areas roughly follow either a r 1 / 4 or an exponential law. 
These properties give support to the young host galaxy interpretation of the extended 
light around QSOs at high-redshift. The rest-frame UV and UV-optical colours are 
inconsistent with the hypothesis of a scattered halo of light from the active nucleus 
by a simple optically-thin scattering process produced by dust or hot electrons. If the 
UV light is indeed stellar, star formation rates of hundreds of solar masses per year 
are implied, an order of magnitude larger than in field galaxies at similar redshifts and 
above. This might indicate that the QSO phenomenon (at least the high-luminosity 
one) is preferentially acompanied by enhanced galactic activity at high-redshifts. 
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1. INTRODUCTION 

Recent deep-imaging studies have discovered a large pop- 
ulation of young galaxies at redshifts z > 0.3, about 3% 
being at z ^ 2.5 for 7? <, 25 (e.g. Steidel et al. 1996a). 
The galaxies have properties very similar to those of present 
day star-forming dwarf galaxies or H II galaxies, but with 
star formation rates SFR<^ 3OM0/yr they are not particu- 
larly luminous objects (Guzman et al. 1997; Phillips et al. 
1997). The field population of high-z galaxies has been as- 
sociated with early counterparts of disks (Wolfe et al. 1995) , 
small spheroids of present-day L J> L* galaxies (Steidel et al. 
1996a) and/or with a collection of star-forming knots being 
merged into galaxies (Lowenthal et al. 1997). On the other 
hand, detailed studies of the intrinsic properties of nearby 
cluster ellipticals indicate that their stellar population was 
formed at z > 2 (Bower, Lucey & Ellis 1992). One long 
standing view is that these were systems formed monolithi- 
cally in an early burst of star formation (e.g. Larson 1974), 
and have been passively evolving since then. Thus, for a first 
rank elliptical we expect SFR~ 10 12 M Q /10 9 yr. Several at- 
tempts to find the early counterparts of these galaxies in 
large areas of the sky have been conducted (see Pritchet 
1994 for a review) but to date just a handful of candi- 
dates that still need spectroscopic confirmation have been 
found (e.g. Pahre et al. 1997). Where are these bright ob- 
jects hiding? In galaxy formation models based on the hier- 
archical merging of dark matter halos the majority of large 
spheroids are formed at lower redshifts (Kauffmann, White 



& Guiderdoni 1993), so that the present-day large systems 
were previously broken-up into small fainter pieces of a yet 
not completely assembled whole. Although they are an at- 
tractive solution, it still needs to be demonstrated whether 
these models can reproduce the intrinsic properties of mas- 
sive ellipticals, such as their large nuclear metal content, the 
thinness of the fundamental plane, the small scatter of the 
colour-luminosity relation, the large colour/metal gradients, 
etc. (but see Kauffmann & Chariot 1997 for first positive re- 
sults on the colour-luminosity relation) . There is also a good 
possibility that some large spheroids may have gone unde- 
tected as the host galaxies of bright QSOs. Two pieces of 
observational evidence support this hypothesis. First, a sig- 
nificant number of nearby bright QSOs are associated with 
luminous ellipticals (McLeod & Rieke 1995; Disney et al. 
1996; Taylor et al. 1996). Secondly, many of the properties 
of high-z QSOs are consistent with their association with 
the cores of early luminous spheroids. Among them, a) the 
high metal content of the Broad Line Region of high-z QSOs 
(Hamman & Ferland 1993); b) the detection of large masses 
of dust in z ~ 4 QSOs (Isaak et al. 1994); and c) the cor- 
respondence of the observed luminosity function of QSOs 
and ellipticals (Terlevich & Boyle 1993 , Haehnelt & Rees 
1993) do suggest that QSO activity might be occurring in 
the central regions of young massive ellipticals. 

All these points are reinforced by the recent discovery 
that about 5 to 20% of the luminosity of both radio- loud and 
radio-quiet QSOs at z ~ 2 arises from extended structures 
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of FWHM « 1-2 arcsec (Lehnert et al. 1992, Aretxaga, 
Boyle & Terlevich 1995, Hutchings 1995). Unraveling the 
nature of the extensions around QSOs at these redshifts is 
of major importance since they could be revealing the hidden 
signature of the elusive early large spheroid population. 

In this paper we focus our attention on the multicolour 
properties of the hosts of three z ~ 2 QSOs previously de- 
tected in i?-band by Aretxaga, Boyle & Terlevich (1995; 
hereinafter ABT95), in an attempt to characterize the emis- 
sion mechanism. In section 2 we describe the dataset. In 
section 3 we describe the method of finding the hosts, and 
their analysis. In section 4 we discuss the nature of the hosts 
on the basis of their colours, luminosities, sizes and profiles. 
In section 5 we summarize our main conclusions. 



The data were reduced with the IRAFj] software pack- 
age. Each QSO frame was first bias subtracted and then 
fiat-fielded using a sky flat-field. The sky-flat field was com- 
puted for each night from the median of the frames of all 
QSO fields, after bright stars, galaxies and QSOs had been 
masked from the data. 

Photometric zero-points for the R and / band data on 
the fields around 1630.5+3749 and PKS 2134+008 were ob- 
tained based on observations made with the 2.5m Issac New- 
ton Telescope (INT*) during service time on 1996 September 
19. Magnitudes were calculated through circular apertures 
which reached the asymptotic value of the circular aperture 
growth, typically about 3 to 4 arcsec in radius. 



2. DATA SET 

The three QSOs studied here were previously detected to be 
extended in the i?-band imaging program of ABT95. They 
comprise the three extended objects in their four target list. 
The QSOs were originally selected from the Veron-Cetty and 
Veron (1994) catalogue with the only condition of belong- 
ing to a narrow redshift-luminosity band (1.8 z <J 2.2, 
M B £ -28 mag for H = 50 Km s" 1 Mpc" 1 and q = 0.5), 
and lying close in projection (20 ^ 9 <, 50 arcsec) to stars of 
similar brightness in order to define a reliable point-spread- 
function (PSF). The QSOs, two radio-quiet and one radio- 
loud, can thus be said to be representative of the high- 
luminosity end of the QSO population at the epoch when 
QSOs were most abundant (e.g. Boyle et al. 1991). 

Our observations (including those already published in 
ABT95) are summarized in Table 1, and detailed in sections 
2.1 and 2.2. 



2.1. Optical Imaging 

The principal optical observations in this study were per- 
formed at the f/11 Auxiliary Port of the 4.2m William Her- 
schel Telescope (WHlQ) in La Palma. The 1024x1024 pixel 
CCD TEK chips in this configuration give a spatial reso- 
lution of 0.105 arcsec/pixel over an unvignetted 1.8arcmin 
diameter field. The CCD was read out in QUICK mode. 

Harris _R-band images for all three QSOs were obtained 
in August 1994 under non-photometric conditions. The re- 
duction and analysis of this data set was already discussed 
in ABT95, and does not differ from that performed in the 
present analysis. In 1995 September 3/4 we acquired Har- 
ris /-band images of the three QSO fields, and a second set 
of _R-band observations for Q 2244—0105. The observations 
were also obtained under non-photometric conditions. 

For each observation, the total integration time of each 
QSO was divided in 600s exposures to avoid saturation of 
the QSOs and stars. Individual frames were offset from one 
another by ~ 10 arcsec to allow discrimination of bad pixels. 

* The William Herschel and the Issac Newton Telescopes are op- 
erated on the island of La Palma by the Royal Greenwich Obser- 
vatory in the Spanish Observatorio del Roque de los Muchachos 
of the Institute de Astrofi'sica de Canarias 



2.2. Near-infrared Imaging 

The near-infrared imaging was conducted at the 3.5m tele- 
scope of the German-Spanish Astronomical Center on Calar 
Alto f| with MAGIC, the 256x256 pixel NICMOS3 CCD ar- 
ray. The f/10 configuration in high-resolution mode gives 
a spatial resolution of 0.32 arcsec/pixel over a square 
82x82 arcsec field. The chip was read out in the standard 
reset. read. read mode. 

On 1995 August 5/6/7 we obtained K' images of the 
three QSO fields under photometric conditions The total 
integration times for the QSO fields were split into 2s units 
to avoid sky saturation, and stacked in 1 min frames. As 
with the optical imaging, individual frames were offset from 
one another by ~ 10 arcsec, making sure that both QSOs 
and nearby PSF stars were included in the field-of-view of 
the detector. 

The data were reduced with the IRAF software package. 
Sky and dark-current were subtracted in each frame using 
the median of six unregistered frames of the same field ob- 
tained just before and after each observation, and scaled up 
to the mean value of source-free patches of the sky. The 
frames were then flat-fielded using a dark-current corrected 
exposure of the uniformly illuminated dome. Each QSO field 
was co-registered and co-added using the centroids of the 
brightest stars and QSOs. Flux calibration was performed 
through a set of standard stars from the list of Elias et al. 
(1982) throughout the nights. Integrated magnitudes for all 
the objects were obtained through circular aperture pho- 
tometry that reached the asymptotic values of the circular 
aperture growth. 



2.3. Near-infrared Spectroscopy 

A 5 min integration spectrum of PKS 2134+008 and a 10 
min integration spectrum of Q 2244—0105 were also ob- 
tained for our program during service time with the CGS4 



t IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Fundation 

t The German-Spanish Astronomical Center on Calar Alto is 
operated by the Max-Plack-Institute for Astronomy, Heidelberg, 
jointly with the Spanish National Comission for Astronomy 
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Table 1. Summary of observations 



Name 



1630.5+3749 



PKS 2134+008* 



Q 2244-0105 



Ml 



-28.3 2.037 



-29.6 1.936 



-28.6 2.040 



Filter 

R 
I 

K' 
R 
I 

K' 
R 
R 
I 

K> 



Exposure 
(seconds) 
10500 
4800 
13320 
5850 
4800 
5220 
5400 
3600 
3600 
7670 



* Mg and z from Veron-Cetty and Veron (1994). 
t Stellar FWHM measured on co-added frame. 

* radio-loud QSO 



seeing* Telescope 
(arcsec) 
0.6 
0.75 
1.0 
0.6 
0.9 
1.1 
0.7 
0.9 
0.7 
1.0 



WHT 4.2m 
WHT 4.2m 
Calar Alto 3.5m 
WHT 4.2m 
WHT 4.2m 
Calar Alto 3.5m 
WHT 4.2m 
WHT 4.2m 
WHT 4.2m 
Calar Alto 3.5m 



Detector 

TEK1 

TEK5 

MAGIC 

TEK1 

TEK5 

MAGIC 

TEK1 

TEK5 

TEK5 

MAGIC 



Date 

1994 Aug 

1995 Sep 
1995 Aug 

1994 Aug 

1995 Sep 
1995 Aug 

1994 Aug 

1995 Sep 
1995 Sep 
1995 Aug 



spectrograph at the 3.8m UK Infrared Telescope (UKIRT|Sj) 
on Mauna Kea. The observations were acquired with spec- 
tral resolution R=550 and spectral coverage 1.45-2.05/x, in 
order to contain the rest-frame [O IIIJA5007 A to Ho region. 
The spectra were flux and wavelength calibrated at the ob- 
servatory following their standard procedures. 

3. DATA ANALYSIS 
3.1. Optical Imaging 

3.1.1. PSF subtraction 

The procedure used to fit and subtract the stellar PSF from 
the images was described in ABT95 in detail. Briefly, it in- 
volved the definition of a two-dimensional PSF from the 
brightest stellar companion to the QSO (marked in Fig. 1 of 
ABT95), and the subtraction of the 2D-profile from the QSO 
and nearby stars, following three recipes: (1) flux-scaled sub- 
traction; (2) subtraction that yields zero flux in the center 
of the QSO; (3) subtraction that yields a smooth flat-top 
profile with no depression in the center. Criteria (2) and (3) 
were evaluated for each frame after smoothing the residu- 
als of the subtraction with a Gaussian filter of a = 1 pixel. 
This procedure provided stability to the solution found from 
frame to frame. The frames were then registered and co- 
added using the centroids of stars and QSOs in the original 
unsubtracted frames. Only non-smoothed frames were used 
in the co-addition. 

Flux-scaled subtraction typically produces solutions 
with negative flux in the center, defined as the centroid of the 
original QSO position, surrounded by a ring of positive flux 
(see Fig. la). For the three QSOs the flux-scaled PSF sub- 
traction is a solution that is inconsistent with QSOs and PSF 
stars having the same spatial profiles. The QSO residuals left 
after the QSO— PSF subtraction depart from the subtraction 
error (a) by over 3cr in both R and I-band images, with the 
only exception being the /-band image of Q 2244—0105 (2a). 

§ The United Kingdom Infrared Telescope is operated by the 
Joint Astronomy Center on behalf of the UK Particle Physics 
and Astronomy Research Council 



As a check, bright nearby PSF-subtracted stars leave residu- 
als around la (see also Figure 2 in ABT95 for a graphic com- 
parison of stellar and QSO residuals in the ii-band frames) . 

Fig. lb and lc show typical radial profiles for the two 
other subtraction criteria. These will be regarded as a lower 
limit and a good approximation to the actual profile of the 
hosts. 

Fig. 2 shows contour plots of the central-zero-flux sub- 
traction solutions. The putative hosts are roundish and show 
no obvious signs of disrupted morphologies for our limit- 
ing PSF-subtraction backgrounds (/j,r « 25.0 to 25.5 mag, 
Hi ~ 23.0 to 23.5 mag at la per pixel, evaluated from the 
sky in the PSF-subtracted area around the QSOs) . 

3.1.2. Accuracy of PSF subtractions 

To assess the accuracy of the PSF subtractions and esti- 
mate the errors on the derived properties of the hosts, we 
have constructed a set of simulated QSOs by superposing 
observed stellar PSFs on compact faint field galaxies and 
also on compact simulated galaxies (FWHM~ 0.8 — 1.2 arc- 
sec). The simulated galaxies were introduced in our images 
with the package ARTDAT of IRAF with the same Poisso- 
nian noise as that of the observations. Nuclear components 
using PSF models derived from field stars and containing 
90, 95, 97, and 99 per cent of the total luminosity of the 
simulated QSOs were superposed on the galaxies, and the 
PSF subtraction was then attempted with the same PSF 
models. The results for two of these galaxies are shown in 
Fig. 3. 

As one would assume intuitively, the simulations show 
that the larger the size and the larger the contribution of the 
host to the total luminosity of the QSO, the easier it is to re- 
cover its true luminosity, radius and profile. The simulations 
also show that flat-top subtractions tend to underestimate 
the total luminosity of the hosts. This effect is actually intro- 
duced by design through the smoothing performed in order 
to evaluate whether a flat regime had been attained. This 
shows that although flat-top PSF subtractions yield better 
estimates for the host properties than central-zero-flux sub- 
tractions, they are still conservative. The biggest departures 
from the true profiles (about a factor of 2) are expected in 
the inner stellar FWHM scale for the range of galaxy com- 
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PSF subtraction methods 
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Figure 1. Radial profiles of the solutions found by performing PSF subtractions on the i?-band frames of 1630.5+3749. Error bars 
represent the scatter around the mean value, and do not include the subtraction error, (a) flux-scaled subtraction; (b) subtraction to 
obtain zero flux in the position of the QSO centroid; (c) subtraction to obtain a smooth flat-top profile 



pactness explored. These discrepancies will be even greater 
for even more compact galaxies. As an example, from the cir- 
cular aperture curves of growth of the most compact galaxies 
in our simulations (FWHM« 0.8 — 1 arcsec), we derive that 
the flux of the hosts can be measured within better than a 
factor of 2 if the host contribution is between 1 and 5 per 
cent of the total QSO luminosity, and increasingly better for 
larger luminosities and sizes. The effective radii of the hosts 
can be derived with relative errors between 20 and 10 per 
cent, for contributions of the galaxy between 3 and 5 per 
cent of the QSO luminosity. 



3.1.3. Properties of the hosts 

Integrated magnitudes, effective radii, and radial profiles of 
the fiat-top PSF subtracted hosts were computed. The re- 
sults are summarized in Table 2 and Fig. 4. 

The magnitudes of the QSO systems (i.e. nuclear + host 
light) are calibrated using our own observations, except for 
the _R-band image of Q 2244—0105, which was calibrated 
using the crude zero-points of ABT95 based on the APM 
Northern Sky Catalogue (Irwin, Maddox & McMahon 1994). 
The difference between these zero-points and our calibration 
for the objects in the other two QSO fields is within the 



error bars of both measurements, but our new calibration is 
typically 0.4 mag fainter. 

Flat-top PSF subtractions yield hosts that comprise 
between 5 and 12 per cent of the total luminosity of the 
QSO+host systems. 

Table 3 contains the colours of the nuclei and host galax- 
ies. The colours of 1630.5+3749 (R — I = 0.8 + 0.1 mag) and 
PKS 2134+008 (R - I = 0.6 ± 0.1) are within the range of 
observed colours of QSOs at that redshift (Veron & Hawkins 
1995). These colours are somewhat bluer than those derived 
for the hosts: R-I ~ 1.1+0.1 mag and R-I ~ 0.8+0.1 mag, 
respectively. However, the errors associated with the PSF 
subtraction are usually larger than the calibration errors, 
and could easily account for this discrepancy between QSO 
and host colours. In order to obtain a more reliable value of 
the colours of the hosts, we have integrated the total lumi- 
nosity enclosed between two concentric circular apertures of 
radius 0.7 arcsec and 3 arcsec. The outer aperture reaches 
the plateau of the aperture growth curves, and the inner 
one neglects the part of the hosts more affected by subtrac- 
tion errors (see Fig. 3). For the host of 1630.5+3749 we find 
R- I = 1.03 and for that of PKS 2134+008 R-I = 1.1, 
redder by 0.24 and 0.51 mag, respectively, than the total 
light of the QSOs. The colours remain similar within rms 
fluctuations of 0.15 and 0.06 mag respectively when the in- 
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PKS2134+008 




Q2244-0105 




Figure 2. R and /-band contour plots of the solutions found after performing a PSF subtraction that yields zero flux in the centroid of 
the QSO. The minimum contour level was set at la of the sky noise after subtraction of the PSF. Further levels indicate lcr increments 
over the minimum level. North is to the left, and East is down. 
© 0000 RAS, MNRAS 000, 000-000 
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Figure 3. Results of flat-top PSF subtractions on simulated QSOs that use galaxies of (a) FWHM= 0.9 and (b) FWHM= 1.2 arcsec, 
with nuclear profiles of FWHM= 0.8 arcsec. The solid lines represent the original profiles of the galaxies before superposing the unresolved 
components. Dashed lines represent the recovered profiles of the galaxies after PSF subtraction. For the simulations presented here, the 
galaxies comprised 3, 5, and 10 per cent of the total luminosity of the simulated QSOs. 



Table 2. Integrated optical magnitudes and scale-lengths of hosts and QSOs. The columns give: name of 
the QSO; filter; m(QSO) magnitude of QSO+host; mi magnitude of the hosts found with central zero-flux 
PSF subtraction; ni2 magnitude of the hosts found with flat-top PSF subtraction ; D e f j effective diameter 
of the hosts found with flat-top PSF subtraction; seeing of the co-added frames, defined as the FWHM of 
the stellar profiles. 



Name 


Filter 


m(QSO) 


mi 


n%2 


D eff 


seeing 






(mag) 


(mag) 


(mag) 


(arcsec) 


(arcsec) 


1630.5+3749 


R 


18.34 + 0.12 


22.0 


21.3 


0.9 


0.6 




I 


17.55 + 0.10 


20.6 


20.2 


1.0 


0.75 


PKS 2134+008 


R 


16.78 + 0.10 


20.8 


20.2 


0.95 


0.6 




I 


16.19 + 0.10 


19.7 


19.4 


1.2 


0.9 


Q 2244-0105 


R 


17.5 + 0.30 t 


m(QSO) +3.4 


m(QSO) +2.3 


0.8 


0.7 




I 




m(QSO) +3.6 


m(QSO) +3.3 


1.0 


0.7 



t Flux calibration based on APM Northern Sky Catalogue. 



ner aperture is set from 0.7 arcsec to 1 arcsec, indicating 
that the outermost regions of the hosts of PKS 2134+008 
and 1630.5+3749 are redder than their nuclei. 

3.2. Infrared Imaging 

A two-dimensional PSF subtraction technique similar to 
that described above for the optical images was also per- 
formed in the infrared frames. The only difference is that 
the individual frames were registered and co-added before 
performing the subtraction. The residuals left in the QSO 
positions were in all cases negligible. 

Radial profile comparison of QSOs and PSF stars in 
the field show that the differences between them are indeed 
small, and located in the outermost regions, where the sky 
subtraction is critical. Fig. 5 shows the radial profiles nor- 
malized to the first computed ring at 0.32 arcsec for QSOs 
and stars. Only 1630.5+3749 has a marginal excess above 
the stellar profile at 2 to 4 arcsec from the center. The ex- 
cess does not disappear even after allowance for shifts in the 
sky level. 

The colours of the QSOs 1630.5+3749 (R - K = 2.5 ± 
0.1 mag), PKS 2134+008 (R - K = 2.1 ± 0.1 mag) and Q 



2244-0105 (R - K » 3.4 ± 0.3 mag) are similar to those 
measured in other QSOs at these redshifts (Hewett, priv. 
communication). If we assume that the excess present in 
1630.5+3749 is real, the surface brightness colour over the 
2-3 arcsec region is R — K — 3.3 + 0.2 mag, redder by about 
0.8 mag than the observed colour of the QSO. 

The surface brightness limits for the K-banA images 
are between 22.1 and 21.8 mag/arcsec 2 at la per pixel. The 
limits for 3<r detections of point-like objects with a typical 
seeing-disk range between 20.8 and 21.1 mag. These limits 
are similar to those obtained in a previous if-band imaging 
program of radio-quiet QSOs at z ~ 2 that yielded non- 
detections (Lowenthal et al. 1995). 

From the K-band limits we can set the colours of the 
hosts of PKS 2134+008 and Q 2244-0105 to be R - K <i 
3.3 mag, consistent with those of the unresolved nuclear 
light. 
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Figure 4. Radial profiles of the hosts derived from flat-top PSF subtractions. The error bars in the inner regions are dominated by the 
error in the subtractions, estimated from the simulations, and in the outer regions by sky subtraction. The /-band profile of Q 2244—0105 
is not flux-calibrated. 



Table 3. Colours of nuclei and hosts. The columns give: name of the QSO; colour; colour of nucleus; total 
colour of host (1) from the flat-top subtractions; colour of the host (2) in the outer 1 to 3 arcsec. 



Name 


Colour 


colour(AGN) 


colour(host) 1 


colour (host) 2 






(mag) 


(mag) 


(mag) 


1630.5+3749 


R- I 


0.8 + 0.1 


1.1 ± 0.1 


1.0 + 0.15 




R-K 


2.5 + 0.1 




3.3 + 0.1 


PKS 2134+008 


R-I 


0.6 + 0.1 


0.8 + 0.1 


1.1 + 0.1 




R-K 


2.1 + 0.1 




£3.3 


Q 2244-0105 


R-I 










R-K 


3.4 + 0.3 




£3.3 



4. DISCUSSION: THE NATURE OF THE 
HOSTS 

The hosts we have detected have luminosities that con- 
tribute between 5 and 12% of the luminosity of the QSOs 
(nucleus+extension) both in R and /-bands (A2300 A and 



A2800 A rest-frame). If the marginal excess over the stellar 
profile of 1630.5+3749 is real, the if-band (A7600 A rest- 
frame) luminosity of the host as derived from the measured 
colours would also contribute to the total A"-band QSO lu- 
minosity about 5%. Our measurements indicate large and 
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Figure 5. iif'-band radial profiles of QSOs and stars in the field, normalized to their central flux. 



luminous extended systems (D e ff ~ 4 Kpc, R w 21 to 
22 mag) . We will consider three alternative explanations for 
their nature: scattered light from the active nucleus, neb- 
ular continuum produced by the gas ionized by the active 
nucleus, and stellar light. 

4.1. Scattering 

The large structures detected in high-redshift radio-galaxies 
(e.g. Chambers et al. 1987) have been attributed to light 
scattered from the active nucleus by a powerful trans- 
verse radio-jet (e.g. Fabian 1989). Indeed, recent spectro- 
polarimetry of two z Ri 1 radio-galaxies that exhibit 'align- 
ment effects' (optical structures oriented in the direction of 
the radio-jet) indicates that around 80% of the total UV con- 
tinuum emission at rest-frame A2800 A is non-stellar scat- 
tered light (Cimatti et al. 1997). Since radio- galaxies and 
radio-loud QSOs could be identical objects viewed from dif- 
ferent angles (Barthel 1989), an important amount of scat- 
tered light may be present around radio-loud nuclei (see Fos- 
bury 1997 for a detailed discussion). 

The host of the core-dominated radio-loud QSO studied 
in this paper, PKS 2134+008, has an R — I colour that 
is 0.52 mag redder than the nucleus. This argues against 
the simple optically-thin scattering case, which would yield 



colours as blue as or bluer than the nucleus itself. The R — K 
limits for this host (<, 3.3 mag) are not deep enough to probe 
the optically-thin scattering case value. Note that four of the 
six lobe-dominated radio-loud QSOs at z ~ 2 — 2.5 studied 
by Lehnert et al. (1992) have hosts with B — K colours 
redder than the nuclei. However the B — K (or R — K) 
analysis discriminates poorly the origin of the light below 
the A4000 A break, since scattered light would contribute 
predominantly in B-band (rest frame A1400 A) and stellar 
light would contribute predominantly in Jf-band (rest-frame 
A7300 A) (e.g. Fosbury 1997). 

The host of our radio-quiet QSO with calibrated pho- 
tometry, 1630.5+3749, exhibits the same general properties, 
with colours 0.25 and 0.8 mag redder than the nucleus in 
R — I and R — K, respectively. These colours are inconsis- 
tent with the optically-thin scattering case as an explanation 
for the UV hosts of our QSOs. An optically-thick medium 
should be invoked in order to produce colours redder than 
those of the scattered source, but then the geometry of the 
scatterers should not be symmetrical since we still see the 
blue colours of the nucleus itself. This could be an alterna- 
tive mechanism for the origin of the extended light in our 
sample, although there is no evidence in general for asym- 
metric scattering in radio-quiet QSOs. 
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4.2. Nebular light 

In ABT95 we argued that nebular continuum produced by 
an extended narrow line region around the active nucleus 
was unlikely to be the origin of the hosts we have detected, 
the reason being that the expected narrow emission lines 
should be very prominent in that case, with peak intensities 
of the Balmer lines more than 3 times larger than those of 
the broad lines. We have since then acquired near-IR spec- 
tra of two of the QSOs studied here (see Fig. 6) which show 
prominent broad Ha lines, but no prominent narrow compo- 
nents. Note that even [O IIIJA5007 A, marked in the figure, 
is marginally detected at best 



4.3. Stellar light 

Stellar light remains still the most plausible interpretation 
for the extended light we have detected. The colours, sizes, 
luminosities and radial profiles are indeed in agreement with 
those expected from young stellar populations: 

a) The colours of the hosts, R — K m 3.3 mag for 
1630.5+3749 and R - K <, 3.3 for PKS 2134+008 and Q 
2244—0105, are not consistent with the colours predicted 
from a simple passively evolved stellar population, usually 
assumed to be characteristic of elliptical galaxies. Bres- 
san, Chiosi & Fagotto (1994), for instance, calculate that 
a passively-evolved stellar population for iz~2 elliptical 
galaxy be R — K ss 5.5 mag and slightly bluer for a mild 
amount of activity. Younger populations are necessary in or- 
der to account for the blue colours observed in our hosts. As 
an example, typical nearby H II galaxies set at z = 2 would 
have colours as blue as R — if ~ 1 — 1.7, as derived from 
their characteristic fiat spectral energy distributions (SED): 
/„ oc v", with <; a <, 0.5. 

b) and c) The luminosities and radii of our hosts lie 
along the luminosity-radius relationship of the local young 
H II galaxies (Telles, Melnick & Terlevich 1997). In Fig. 7 
we converted the UV luminosities of the hosts (observed R- 
band) to rest-frame B-band using the SED of local H II 
galaxies. This is equivalent to converting the £?-band lumi- 
nosities of the local sample of H II galaxies to rest-frame 
A2300 A, and then comparing them with those of the hosts. 
Note that there is at least one local H II galaxy that is as 
big and luminous as our hosts. 

d) The radial profiles of the i?-band hosts, derived from 
the flat-top solutions, fall approximately as r 1 ^ 4 -laws or ex- 
ponential profiles for radii r J> 0.6 arcsec (Fig. 8). Profiles 
derived for radii smaller than the FWHM of the observations 
are usually unreliably recovered by flat-top subtractions, as 
shown by our numerical simulations of galaxy+PSF (see sec- 
tion 3.1.2). Distinguishing between different morphologies is 
a hopeless task at these redshifts since the discrimination 
between disks and bulges is most efficiently carried out in 
the the inner part of the profiles, where our data are unreli- 
able. The point we want to illustrate with Fig. 8 is that the 
morphology of the hosts is not definitely different from that 
of smooth local galaxies, at least in the outer regions. 



5. CONCLUSIONS 

We have shown that the hosts of z ~ 2 high-luminosity 
radio-quiet and radio-loud QSOs are large and luminous, 
rivaling the most luminous nearby star forming galaxies. 
Colours, sizes, luminosities and radial profiles are compati- 
ble with the interpretation of the extended light as a young 
galaxy. If all the UV luminosity is coming from a stellar 
population, the implied star formation rates (in the contin- 
uous star formation case) are SFRJ> 100 — 200 Me/yr on 
scales of D e ff « 4 Kpc. These values are about an order 
of magnitude higher than those derived for field galaxies at 
similar redshifts found by Lyman Break searches in the Hub- 
ble Deep Field (Steidel et al. 1996b, Lowenthal et al. 1997). 
At z ~ 2, an unevolved L* galaxy with SED typical of a star- 
forming galaxy would appear to be about 3 mag fainter than 
the hosts we have detected. The density of high-luminosity 
QSOs (Mb —28 mag) at redshifts between 1 <, z <, 3 
(Boyle et al. 1991) is about 10 Gpc -3 . This means that the 
density of large luminous galaxies like the ones considered 
here should be about 8 Gpc -3 at 1 < z < 3, taking into 
account the non-detection case of ABT95 and the relative 
under- abundance of radio-loud QSOs, which form less than 
10% of the QSO population. 

While these results are intriguing, it is still possible that 
other sources such as nebular light and scattering may con- 
tribute to the observed extensions. However, the absence 
of prominent narrow emission lines in the spectra of these 
QSOs argues against a major contribution by the nebular 
continuum emitted by an extended narrow line region. The 
optically-thin scattering hypothesis is furthermore ruled out 
as a major mechanism by redder UV and UV-optical colours 
than those of the active nuclei. The optically-thick scatter- 
ing case is still an open possibility, especially in radio-loud 
hosts, where a role for scattering has already been invoked. 
High spatial resolution spectroscopic observations of the ex- 
tended light is probably the simplest and most direct way 
to assess the contribution of all those different mechanisms 
into the UV luminosities we have measured. This is an im- 
portant project in order to reliably prove and measure the 
contribution of massive star-formation in QSO hosts. 

The study of the hosts of less luminous QSOs, espe- 
cially those which populate the break of the QSO luminosity 
function, is also of major importance, in order to assess the 
association of large luminous galaxies and the majority of 
QSOs at the QSO epoch. 
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